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We report on a atomistic theory of electronic structure and optical properties of a single InAs 
quantum dot grown on InP patterned substrate. The spatial positioning of individual dots using 
InP nano-templates results in a quantum dot embedded in InP pyramid. The strain distribution 
of a quantum dot in InP pyramid is calculated using the continuum elasticity theory. The electron 
and valence hole single-particle states are calculated using atomistic effective-bond-orbital model 
with second nearest-neighbor interactions, coupled to strain via Bir-Pikus Hamiltonian. The optical 
properties are determined by solving many-exciton Hamiltonian for interacting electron and hole 
complexes using the configuration-interaction method. The effect of positioning of quantum dots 
using nanotemplate on their optical spectra is determined by a comparison with dots on unpatterned 
substrates, and with experimental results. The possibility of tuning the quantum dot properties with 
varying the nano-template is explored. 
PACS numbers: 73.21.La, 78.67.Hc, 71.15.-m 



I. INTRODUCTION 

There is currently significant interest in semiconduc- 
tor self-assembled quantum dots (SAQDs) 0, 0, 13 due 
to their excellent electronic and optical properties. The 
quality of optical properties is primarily due to the very 
clean self-organized Stransky-Krastanow growth process 
during molecular beam epitaxy. The downside of this 
process is the random spatial distribution of quantum 
dots. In order to combine single quantum dots with cav- 
ities, gates, or magnetic ions for their increased function- 
ality one needs to position quantum dots by growth on 
patterned substrates . 

The patterned substrate growth is both a challenge and 
opportunity. It is a challenge to assure that patterning of 
the substrate does not destroy the high quality of single 
quantum dots, and it is an opportunity to use patterning 
to control their electronic properties. Recently Williams 
et al 0, ||| reported high quality emission spectra close 
to 1.55 fim of single InAs quantum dots grown on InP 
pyramidal nanotemplate. This opens up the possibility 
of integration of single quantum dots into optical cavi- 
ties 0, and hence reliable fabrication of single photon 
sources @, suitable for long distance transmission in 
optical fiber. The theoretical understanding of the ef- 
fect of patterning on optical spectra is unknown, and we 
present here atomistic theory of electronic and optical 
properties of single InAs quantum dot grown on pyrami- 
dal InP nanotcmplates. 

Atomistic calculation of the electronic properties of 
InAs dots on InP presents a challenge: due to small lat- 
tice mismatch between InAs and InP (3%), the InAs 
dots are significantly larger then InAs dots on GaAs. 
The larger dot size would suggest the applicability of the 
k • p method, unfortunately it suffers from the unphys- 
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FIG. 1: Single InAs quantum dot on a pyramidal InP nan- 
otemplate (a) before overgrowth, enlarged view of top portion 
of the structure and (b) after overgrowth. 



ical level crossing as pointed out by Holma et al jlfij . 
Hence atomistic method, which gives proper bulk band 
dispersion, is needed for InAs/InP dots. In this work, 
the electron and valence hole single particle states are 
calculated using atomistic effective-bond-orbital model 
(EBOM) 0, Il2f with second nearest-neighbor interac- 
tions, coupled to separately calculated strain distribution 
via Bir-Pikus Hamiltonian |13|. The optical properties 
of InAs dots embedded in InP pyramids are calculated 
by solving the many-exciton Hamiltonian for interact- 
ing electron and hole complexes using the configuration- 
interaction method. The effect of positioning of quan- 
tum dots using nanotemplate on their optical spectra is 
determined by a comparison with dots on unpatterned 
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substrates, and with experimental results. The possibil- 
ity of tuning the optical properties with the shape of the 
nanotemplate is demonstrated. 



II. MODEL OF INAS QUANTUM DOT ON 
NANOTEMPLATE 

In the work of Williams et al 0, 0] the positioning of 
InAs dots is achieved by first the growth of InP nano- 
template with controlled size and shape, followed by the 
growth of a single InAs dot, as shown in Fig. ^ The 
quantum dot is covered by continuing the growth of InP, 
with the final result of a single InAs dot embedded in 
an InP pyramid. There is experimental evidence p| that 
the shape and size of the template determines the ge- 
ometry of the dot. Hence for a rectangular template, 
including facets associated with [100] and [111] crystal- 
lographic axis shown in Fig. ^ we model the dot as a 
rectangular InAs box, with truncations around its four 
corners. 

In this calculation, we chose the size of the template as 
D=A6 nm, the lateral size of the dot as <i=36 nm, with 
height of h=2.35 nm [l4|. After overgrowth, the dot ends 
up embedded 18 nm below the top of the pyramid, as 
shown in the lower panel of Fig. ^ The size and the 
shape of the nanotemplate (denoted by W in Fig. [IJ can 
be varied. The pyramid fabricated in Refs. 0, U had 
W = 400 nm. In this exploratory example, we chose a 
smaller size of the pyramidal template, W = 206 nm with 
height h — 103 nm to reduce significant computational 
effort. 



III. ATOMISTIC TIGHT-BINDING APPROACH 
TO ELECTRONIC STRUCTURE 

Electronic structure of InAs/InP dots has been calcu- 
lated in the past using the quasicontinuous eight-band 
k • p method [lq. For atomistic calculations, the em- 
pirical tight-binding method (ETBM) ^(| is a natural 
choice. However, considering that InAs/InP dots are gen- 
erally much larger than InAs/GaAs dots due to smaller 
lattice mismatch, we chose EBOM due to its less re- 
quirement of computational resources and a compatibil- 
ity with the k * p method. EBOM is a sp 3 tight-binding 
method in which the full symmetry of zinc-blende lattice 
is reduced to that of a fee lattice. The one electron tight- 
binding Hamiltonian describes an electron hopping from 
atomic orbital a at position R to orbital a' at position 
R': 

H tb = ff(R,a;R',a') c R,a c R',«', (1) 

where cj^ q (cr/ iQ ') are creation ( annihilation ) operators. 
The hopping matrix elements and site energies for s and 
p orbitals, extended to include second nearest-neighbor 
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FIG. 2: Band structure of InP described by the eight-band 
k ■ p method (dotted lines) and its modified version (dash 
lines, see text), compared with that by the tight-binding-like 
EBOM (solid lines). 



interactions, are given in real space by 

H(Rp,R'p) = fe-R/.r^^+^a-^)] 

+ E xx 5b.,w, 
H(Rp,R' P ') = K^°t p v<5 r _ r / )T) 

H(TLs,R'p) = E^Vr-RV, (2) 

where r and a give positions of the nearest and second- 
nearest neighbors, respectively, 

r = |[(±1,±1,0),(±1,0,±1),(0,±1,±1)], 

a = a[(±l,0,0),(0,±l,0),(0,0,±l)], (3) 

a is the lattice constant. 

There are ten fitting parameters, E™°, E]f , E™°, 
E^, E™, EH , Elf, Elf in EBOM. The 

two-center approximation |17| introduces an additional 
constraint, E xx ° - E xy a — E° xx x . With second nearest- 
neighbor interactions included 0, effective masses of 
electrons and holes at the T point [ToL IT^ ] , conduction 
and valence band edges at both the T and X points [2(j 
can be exactly fitted within EBOM. With increasing k, 
it can reproduce correct band dispersion along the V — X 
direction as shown in Fig. [21 

Compared with the eight-band k ■ p method, one of ad- 
vantages of the EBOM is the absence of spurious cross- 
ing between valence bands (see Fig. (2J. To alleviate 
this deficiency in the k • p method, terms proportional 
to — (k x + ky + kj) are proposed [T(| to be added to 
heavy-hole and light-hole bands to eliminate the result- 
ing spurious solutions (see dash line in Fig. |2J) . Although 
this modification does not alter effective masses at the T 
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point, it is nevertheless seen to give unphysical dispersion 
of valence bands along the T — X and T — L directions, 
compared with that given by the EBOM. 

When EBOM is applied to two different materials, a 
band offset between the unstrained island and matrix 
material is needed to obtain all the fitting parameters. 
The effect of strain is incorporated into EBOM through 
the deformation potential Bir-Pikus's theory |l2|,[l3| and 
through the piezoelectric effect. An advantage of the 
EBOM over ETBM is its lack of ambiguity in fitting the 
deformation potential parameters |22|. 



IV. STRAIN DISTRIBUTION 

The strain in the vicinity of quantum dot is deter- 
mined by the position of the quantum dot in the pyramid. 
Hence the domain of strain calculation should include 
the entire pyramid, with characteristic sizes on the mi- 
cron scale. While the full atomistic and finite element 
calculation for the quantum dot and the pyramid will be 
implemented in the future, we report here calculations 
based on continuum elasticity theory j2^|. We enclose 
the template in a rectangular computational box. A fixed 
boundary condition is applied to the base of the template 
while free-standing boundary conditions are implemented 
to all the other exposed facets. This is achieved by in- 
troducing dummy sites around the template where in- 
finitesimal elastic constants are assigned, therefore, the 
force applied to the template facets can be neglected. 
The quality of results of continuum elasticity theory cal- 
culations are successfully verified by comparison with the 
results of atomistic valence- force-field calculations :24] for 
smaller structures. 

The strain tensor e is obtained by minimizing the fol- 
lowing elastic energy functional 
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FIG. 3: Plot of hydrostatic and biaxial strain in the InAs/InP 
quantum dot along the [001] direction towards the apex (a) 
and along [100] direction (b). For comparison, strain profiles 
are shown in dash lines for a InAs dot of the same geometry 
in an unpatterned InP substrate. 
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where Cn, C44, and C12 are elastic constants of quantum 
dot material and of InP pyramid material, the values of 
which are taken from Ref. jlOj . 

Figure 13 shows two typical strain components, hydro- 
static and biaxial strain, through the central axis of the 
dot ([001] direction) and through the center of the dot 
along the base ([100] direction). The hydrostatic and 
biaxial strain are defined [19j by 
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The former mainly affects the conduction bands while 
the latter only affects the valence bands, inducing the 
splitting between the heavy-hole and light-hole band. 



The same strain profiles are also shown for InAs dot 
in an unpatterned InP substrate. In this comparison we 
assumed the same shape and size of the quantum dot, 
excluding the possibility that the unpatterned dots are 
disk-like or lens-shaped. Hereafter, dot on unpatterned 
substrate is referred to as dot 1 and the dot on patterned 
substrate is referred to as dot 2. 

The comparison between them shows the effect of pat- 
terned substrate on the strain distribution. Along the 
growth direction [see Fig. Ota)], both the hydrostatic and 
biaxial strain are found marginally larger in dot 2 than in 
dot 1, which would lead to a small blue shift of emission 
energy. Along the direction towards the apex, due to the 
physical limitation of the patterned substrate, obvious 
discrepancy can be found between the strain profiles in 
the two dots that are in different environment. Much less 
difference is found along the opposite direction where the 
patterned substrate becomes similar to bulk InP. Along 
the [100] direction, the strain profiles are found similar 
for the two dots except for their slightly different peak 
values. 
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FIG. 4: Energy levels of the single In As dot on unpatterned 
InP substrate (a) and on patterned substrate (b) calculated 
by the EBOM. It is noted that the energy scale in conduction 
band is different from that in valence bands. 



According to the deformation potential theory ^| , the 
local band edges are proportional to the hydrostatic and 
biaxial strain. If we choose the valence band offset be- 
tween unstrained InAs and InP as 400 meV the 
depth of confinement along the growth direction is 420 
meV and 530 meV for electrons and holes, res pec tively, 
which is in agreement with previous calculation [lCj . The 
heavy-hole band is the top-most valence band, separated 
from the light-hole band by 150 meV due to the biax- 
ial strain. Although the band offset between unstrained 
InAs and InP is not well known, it has been shown that 
the actual electronic structure are not sensitive to the 
value chosen in the calculation [22|. 



V. ELECTRONIC STRUCTURE 

The computational box for electronic calculation is 
chosen to be much smaller than that for the strain calcu- 
lation as we are interested only in the states confined in 
the dot. It contains the whole dot and some sites outside 
the patterned substrate where much larger site energies 
are assigned. Its dimensions are 47.5 nm x 47.5 nm x 
14.1 nm. In EBOM, each unit cell consists of four ef- 
fective atoms, and each atom has eight spin-orbitals. In 
total, it gives rise to a large sparse matrix of dimension 
exceeding five million. The electron states are obtained 
by solving the sparse matrix using the Lanczos algorithm 
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Figure^] (a) and (b) shows the calculated energy levels 
of electrons and holes by EBOM for dot 1 and 2, respec- 
tively. The calculated energy shells are characteristic to 
the square geometry of the dot, with a single s-level, al- 
most degenerate p-shell, and a d-shell consiting of three 
states: a single level followed by two degenerate levels. 




FIG. 5: Probability density of confined electron (upper row) 
and hole (lower row) states in the dot. 



The order and degeneracies (1,2) of levels in the d-shell 
differentiates the square dot from a lens-shaped (3) or 
disk-like quantum dot (2,1) studied previously and 
could be used to infer the shape from the optical spec- 
tra. 

As can be seen, strain profile specific to nanotemplate 
(see Fig. |3J) induces shifts of energy levels for the dot on 
patterned substrate. The shift in the conduction band, 
which is about 11 meV, is more homogeneous than that 
in the valence bands (about 3 meV) as the latter is also 
under the influence of other strain components. 

The probability density of confined states is shown in 
Fig. 03 The hole states are found generally more confined 
than their electron counterparts. Also note that the elec- 
tron and hole states of p-shell are oriented in the oppo- 
site directions. This is a result of piezoelectric potential 
which attracts electrons while simultaneously repelling 
holes. 

We can label the p-like states according to the exten- 
sion direction of their wave functions. In the conduction 
band, the two p-like states are p e + and p e _ , localized along 
[110] and [110] directions, with energies E v ^ < E p « . In 

the valence bands, the two p-like states, and have 
a different order, i.e. E p h < E p h [25J. 

This inverted order of energy of states with similar 
symmetry is due to the enhanced piezoelectric effect in 
the dot on patterned substrate. The piezoelectric poten- 
tial is calculated from the piezoelectric charge density as 
given by 0] 



p P (r) = -2e 
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de yz (r) de xz (v) de xy (r) 



dx 



dy 



(6) 



where en is the piezoelectric constant. Unlike 
InAs/GaAs dots where ei4 in GaAs is four times larger 
than that in InAs, ei4 in InP is smaller than in InAs. 

The left two plots in Figure show the distribution of 
piezoelectric charge density for the two dots in different 
environment. The piezoelectric charge in dot 1 is seen to 
have a nearly symmetric distribution, i.e., a part of pos- 
itive charge appears in pair with another one of negative 
charge. The two parts are almost symmetric and hence 
compensate their effect on electron and hole states, which 
can be seen in the corresponding piezoelectric potential 
shown in the upper-right panel. The potential is found 
to be well localized around the four corners of the dot. 
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FIG. 6: Effect of nanotemplate : isosurface plots of piezoelec- 
tric charge density in the dot on patterned substrate (lower- 
left plot) and the one on unpatterned substrate (upper-left 
plot). Bright (dark) grey area has positive (negative) charge. 
The right plots show the corresponding piezoelectric poten- 
tials as seen by an electron along the quantum-dot layer which 
is about 0.6 nm above the bottom of dots. The peak values 
of both potentials are about ±8 meV. 



Because the positive charge is paired with the negative 
charge, the potential is very similar to that of dipoles, 
and is almost zero inside the dot. This leads to only 0.1 
meV splitting between the p e + and p e _ states. 

In dot 2, the piezoelectric charge density is seen to be 
localized below the dot as well as along the facets of the 
substrate. This asymmetric distribution of the charge 
density induces a very different piezoelectric potential. 
As seen in the lower-right part of Fig. H3 the potential 
inside the dot has significant value. This leads to a split- 
ting of about 0.5 meV between p e + and pi, and inverts 
the order of states in valence bands. 

Figure [7| shows the intensity of interband transitions 
(joined optical density of states) between the calcu- 
lated hole and electron states. Here, we only plot the 
transitions polarized within the plane perpendicular to 
the growth direction, which corresponds to those from 
the heavy-hole part of the hole states. After the ground 
state transition, the next two strong transitions are those 
between the p-like hole and electron states. The transi- 
tions between electron and hole states that have different 



and p] 



PI 



are seen to have 



symmetry, i.e., p_ 
lower intensity. 

A further calculation shows that the tight-binding re- 
sults can be fairly well approximated by a single-band 
effective-mass approach [21| using the following effec- 
tive mass parameters, m* = 0.065, to° 01 = 0.32, and 
= 0.17. The comparison of these parameters with 
those of the bulk materials exhibits the effect of strain 
on the effective mass renormalization in quantum dots. 
For example, the electron effective mass is enhanced from 
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FIG. 7: Intensity of interband transitions between hole and 
electron states. 



0.024 in bulk InAs to 0.065, which is close to the effec- 
tive mass in bulk InP (0.077). The anisotropy in the hole 
effective mass tensor is reversed, the holes become much 
lighter (comparing with 0.639 in bulk InAs and 0.933 in 
bulk InP) within the plane perpendicular to the growth 
direction. This results in comparable energy separations 
in conduction (14.8 meV) and valence (8.0 meV) bands. 



VI. EMISSION SPECTRA OF DOTS ON 
PATTERNED AND UNPATTERNED 
SUBSTRATES 

We now turn to investigate the effect of nano-template 
growth on electronic states as could be observed from 
emission spectra as a function of the population of the 
quantum dot with electrons and holes. The interaction 
of electrons and holes significantly changes the emission 
spectra from those expected from the joint optical density 
of states. The electrons and hole interaction is described 
via the following Hamiltonian, 
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where Ef and E^ are the energy levels shown in Fig. 
0] and V ee , V , and V he are the Coulomb matrix ele- 
ments which can be calculated using the single-particle 
states as shown in Fig. In the calculation reported be- 
low, the multi-exciton configurations arc built from the 
first twelve electron and twelve hole single-particle states 
(with quasi-spin). These configurations describe the low- 
est s, p, and d shells. The Hamiltonian is then expanded 
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FIG. 8: Emission spectrum from individual multi-exciton 
complex obtained by solving the many-body Hamiltonian 
using the configuration-interaction method, shown as peaks 
filled in dark color. For comparison, the results for the dot 
on an unpatterned substrate are shown as unfilled peaks. 
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on these configurations and solved by the configuration- 
interaction method. For more detailed description of this 
calculation we refer to our previous work 21]. 

Figure |S| shows the emission spectra for dot 1 and 2 
as a function of photon energy and increasing number of 
excitons . The emission intensity is calculated for all the 
possible transitions between the n-exciton complex and 
(n-l)-exciton complex at a temperature of 4.2 K. The 
comparison of the spectra between the two dots should 
show the effect of nano-template, in particular the effect 
of inverted order of states in the conduction and valence 
band of the dot on patterned substrate. The inverted 
order of p-like states in valence bands is due to different 
strain distribution, especially in the region close to the 
corners of the dot structure, and hence different piezo- 
electric potential. However, we find very little difference 
in the spectra between the two dots except for an over- 
all blue shift which reflects the difference found in the 
corresponding single-particle energy spectra (see Fig. |4"|. 

Let us first look at tri-exciton which is the first non- 
trivial case. Fig. illustrates four configurations of the 
p-shell states for a tri-exciton with definite spin projec- 
tion in its upper panel. The ground state of the non- 
interacting tri-exciton complex is given by the configu- 
ration \c) = |p_)|f>i) as it has the lowest kinetic en- 
ergy. The other two configurations \a) = \p_)\pt.) and 
\b) = \p+)\p+) have higher kinetic energies due to the 
splittings between the and p_ states. Hence the 



jb2 



la) + lb) 

FIG. 9: Upper panel: Four possible configurations of the p- 
shell states whose probability densities are plotted in the left, 
Lower panel: Schematic energy levels. 



ground state of the non-interacting tri-exciton would be 
dominated by |c), and is not be optically active. The en- 
ergies of tri-exciton configurations and a schematic rep- 
resentation of their oscillator strength (solid line - high, 
dashed line - low ) are shown in the lower panel of Fig. 

M 

The non-interacting tri-exciton ground state is the 
dark configuration |c). However, the order of energy lev- 
els changes once electron-hole attraction (Hartree-Fock 
energy) in each configuration is included. As the spatial 
overlap between p^_ and pi is much smaller (see Fig. 
than either that between p_ and pi, or that between 
p\ and p+, the corresponding Coulomb interaction is 
smaller by 2.3 meV compared with those of the other 
two pairs. Taking this into account, configuration |c) has 
higher Hartree-Fock energy than \a) and \b) by about 1.8 
meV. Configuration \d) has an even higher energy than 
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c) because of the splitting between the p+ and 
states. This arrangement of HF energy levels is shown 
in lower-central panel of Fig. |5J Once correlation are in- 
cluded ( lower-right panel of Fig. [5]), the ground state of 
tri-exciton primarily consists of a symmetric combination 
of configurations \a) and \b) j2||, and becomes optically 
active. 

When the states beyond the second-shell are included, 
the ground state of four-exciton complex (A4) is domi- 
nated by the singlet-singlet (SS) configuration |26J where 
the two electrons (holes) occupying the two p-like states 
have anti-parallel spins. However, a state dominated by 
a triplet-triplet (TT) configuration is lying very close to 
the ground state. Hence, there are two emission peaks 
found in the p-shell emission spectrum of 4X. Depend- 
ing on the temperature and spin relaxation, the emission 
from this excited state may become stronger. The ground 
state of five-exciton complex is well isolated from excited 
states. The two emission peaks found in the p-shell are 
the result of recombinations to the ground and first ex- 
cited 4X states. The emission spectrum from s and p 
shells for dot on patterned substrate is very similar to 
the dot on unpatterned substrate. Hence patterning did 
not deteriorate optical properties. 



A. Power-dependent emission 
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FIG. 10: The calculated photo-luminescence spectra at vari- 
ous excitation power P (in arbitrary units). The major indi- 
vidual emission peaks are labeled according to their origin. 



In the experiment |5j, the emission spectra are mea- 
sured as a function of excitation power and include the 
contribution from all the multi-exciton complexes. In our 
calculation, we solve the rate equations |27| for given ex- 
citation power to obtain the distribution probability of 
each multi-exciton complex. The overall emission spec- 
tra are computed by the summation of the contribution 
from individual multi-exciton complex weighted by the 
calculated probabilities. 

Figure ITUI shows the dependence of photoluminescence 
spectrum on the excitation power controlled by dimen- 
sionless parameter P. The first emission line appearing 
in the s shell is from the single exciton and the second 
line is from a bi-exciton. The calculated single exciton 
emission line is about 15 meV below what the experi- 
ment measured. The calculated bi-exciton binding en- 
ergy (measured by the separation between the X and X2 
lines) is 1.4 meV, which is larger than what the experi- 
ment measured (0.9 meV). 

In the s shell between the X and X2 lines lies an emis- 
sion line from the tri-exciton. As the excitation power 
increases, the emission from the p shell is seen to have 
larger intensity than the s shell. In the p shell, the major 
emission peaks from A3 (A5) and A4 (Ag) overlap each 
other due to their very similar emission energies (see also 
Fig. |SJ. In addition, emission from higher- lying states 
beyond the p-shell can also be seen at higher energies, 
which is characteristic to that from a single 2s level of 
the d-shell, populated by one and two electron-hole pairs, 
much like an s-shell. 



The energy separation between the emission lines of 
the single exciton and tri-exciton defines the separation 
of the s and p shell. Compared with the experimental 
value ,5j, which is about 18.0 meV, our calculation gives 
17.9 meV. This shell separation is largely determined by 
the lateral size of the dot and not sensitive to the height 
of the dot. 



VII. TUNING OPTICAL PROPERTIES WITH 
NANOTEMPLATE 

We now turn to the possibility of tuning the shape of a 
quantum dot with the shape of the template. We assume 
that the change in the template leads to the controlled 
elongation of the quantum dot. The degree of elongation 
is defined as rj = (d — d')/d where d and d' denote the 
dimensions of the rectangular base of the elongated dot. 
We keep d as a constant (36 nm) and change d' to obtain 
different elongated geometries. 

To explain the effect of elongation on the emission 
spectrum we discuss here in detail the emission from the 
five-exciton complex. Even without elongation the emis- 
sion from the 5X complex results in two emission lines 
in the p-shell, as shown in Fig|S] These two lines orig- 
inate from two different 4X states: singlet-singlet (SS) 
and triplet-triplet (TT) 4X states The splitting be- 
tween the SS and TT states is mainly determined by the 
competition between exchange energy of triplets and cor- 
relations in singlet-singlet configurations. When there is 



8 



X 5^ X 4 




1 


| Tl = 10% 


... II II 11 = 6.6% 




I 11 n = 3.3% 




ll 


ss 


ti = 





0.855 0.86 0.865 

Energy (eV) 



FIG. 11: Emission spectra of the five-exciton complex in 
quantum dots with increasing elongation 

elongation in the geometry of the dot, the almost degen- 
erate single-particle states in the p-shell have an extra 
splitting proportional to the elongation. This splitting is 
present in the splitting of the SS and TT states of 4X 
complex, which can be seen in the emission spectrum of 
5X. This is illustrated in Fig. ^2 which shows the de- 
pendence of emission spectra on the elongation of the 
structure. 

The splitting between the two emission peaks of 5X is 
seen to increase with the degree of elongation, which can 
be well explained by the increasing splitting between the 
single-particle states in the p-shell. The emission peaks 
of all the multi-exciton complexes are found to have blue 
shifts as the degree of elongation increases, which is due 
to the shrinking of the volume of the quantum dot. 

We summarize the effect of the elongation on the emis- 
sion spectra as a function of excitation power in Fig. IT51 
which shows the emission spectra for an elongated dot 
with -q = 6.6%. While the elongation is seen not to affect 
the emission in the s-shell, the splittings it induced in the 
p-shell can be clearly identified. The order of individual 
emission peaks in the p-shell is also found sensitive to the 
elongation. 

VIII. CONCLUSIONS 



terned substrate, we have shown that the patterned sub- 
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FIG. 12: The photoluminescence spectra at various excita- 
tion power P for an elongated dot (77 = 6.6%). 



strate breaks the reflection symmetry of the structure, 
enhances the piezoelectric effect, and leads to the in- 
verted order of the p-shell states. We have shown that 
Coulomb interactions restore the order of electron hole 
states and result in emission spectra unaffected by the 
patterning process. We have identified features in the 
emission spectra which can be related to the shape of 
the quantum dots, in particular the characteristic emis- 
sion pattern of the d-shell. We have explored the possibil- 
ity of tuning the emission spectra by changing the shape 
of the quantum dot in a controlled way using the nano- 
template. The theoretical results support the notion that 
the nanotemplate growth of single self-assembled quan- 
tum dots is a promising route for the control of the po- 
sition of self-assembled dots, and hence enables their in- 
creased functionality when combined with cavities, mag- 
netic ions, doping and gates. In turn, it is hoped that 
better experimental control will allow a systematic com- 
parison of microscopic theories with experiment. 



In conclusion, we have presented a theoretical study of 
electronic and optical properties of a single InAs dot on 
InP patterned substrate. We have calculated the strain 
distribution in nanotemplate with appropriate bound- 
ary conditions. By applying the atomistic tight-binding 
effective-bond-orbital method, we have obtained the elec- 
tronic structure of the dot on patterned substrate. By 
comparing the results with those for the dot on unpat- 
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